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Abstract 



We study the average case performance of multi-task Gaussian process (GP) re- 
gression as captured in the learning curve, i.e. the average Bayes error for a chosen 
task versus the total number of examples n for all tasks. For GP covariances that 
are the product of an input-dependent covariance function and a free-form inter- 
task covariance matrix, we show that accurate approximations for the learning 
curve can be obtained for an arbitrary number of tasks T. We use these to study 
the asymptotic learning behaviour for large n. Surprisingly, multi-task learning 
can be asymptotically essentially useless, in the sense that examples from other 
tasks help only when the degree of inter-task correlation, p, is near its maximal 
value p = 1. This effect is most extreme for learning of smooth target functions 
as described by e.g. squared exponential kernels. We also demonstrate that when 
learning many tasks, the learning curves separate into an initial phase, where the 
Bayes error on each task is reduced down to a plateau value by "collective learn- 
ing" even though most tasks have not seen examples, and a final decay that occurs 
once the number of examples is proportional to the number of tasks. 



1 Introduction and motivation 

Gaussian processes (GPs) [ 1 1 have been popular in the NIPS community for a number of years 
now, as one of the key non-parametric Bayesian inference approaches. In the simplest case one can 
use a GP prior when learning a function from data. In line with growing interest in multi-task or 
transfer learning, where relatedness between tasks is used to aid learning of the individual tasks (see 
e.g. 12] 0), GPs have increasingly also been used in a multi-task setting. A number of different 
choices of covariance functions have been proposed J4j [5] [6l [7j [8) . These differ e.g. in assumptions 
on whether the functions to be learned are related to a smaller number of latent functions or have 
free-form inter-task correlations; for a recent review see (9). 

Given this interest in multi-task GPs, one would like to quantify the benefits that they bring compared 
to single-task learning. PAC-style bounds for classification [2 3 1 ) in more general multi-task sce- 
narios exist, but there has been little work on average case analysis. The basic question in this setting 
is: how does the Bayes error on a given task depend on the number of training examples for all tasks, 
when averaged over all data sets of the given size. For a single regression task, this learning curve 
has become relatively well understood since the late 1990s, with a number of bounds and approxi- 
mations available ITT] [T2 [131 [14] [T31 [161 [17] [13 as well as some exact predictions EO). Already 
two-task GP regression is much more difficult to analyse, and progress was made only very recently 
at NIPS 2009 ETI . where upper and lower bounds for learning curves were derived. The tightest of 
these bounds, however, either required evaluation by Monte Carlo sampling, or assumed knowledge 
of the corresponding single-task learning curves. Here our aim is to obtain accurate learning curve 
approximations that apply to an arbitrary number T of tasks, and that can be evaluated explicitly 
without recourse to sampling. 
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We begin (Sec. [2]i by expressing the Bayes error for any single task in a multi-task GP regression 
problem in a convenient feature space form, where individual training examples enter additively. 
This requires the introduction of a non-trivial tensor structure combining feature space components 
and tasks. Considering the change in error when adding an example for some task leads to partial 
differential equations linking the Bayes errors for all tasks. Solving these using the method of 
characteristics then gives, as our primary result, the desired learning curve approximation (Sec. [3]). In 
Sec.fflwe discuss some of its predictions. The approximation correctly delineates the limits of pure 
transfer learning, when all examples are from tasks other than the one of interest. Next we compare 
with numerical simulations for some two-task scenarios, finding good qualitative agreement. These 
results also highlight a surprising feature, namely that asymptotically the relatedness between tasks 
can become much less useful. We analyse this effect in some detail, showing that it is most extreme 
for learning of smooth functions. Finally we discuss the case of many tasks, where there is an 
unexpected separation of the learning curves into a fast initial error decay arising from "collective 
learning", and a much slower final part where tasks are learned almost independently. 

2 GP regression and Bayes error 

We consider GP regression for T functions f T (x), r = 1,2, ... ,T. These functions have to be 
learned from n training examples (xg,Tg,yg), £ = 1, ...,n. Here xg is the training input, Tg g 
{1, . . . , T} denotes which task the example relates to, and yi is the corresponding training output. 
We assume that the latter is given by the target function value f Tt {xt) corrupted by i.i.d. additive 
Gaussian noise with zero mean and variance . This setup allows the noise level a\ to depend on 
the task. 

In GP regression the prior over the functions f T (x) is a Gaussian process. This means that for any 
set of inputs xg and task labels ti, the function values {f T{ (xi)} have a joint Gaussian distribution. 
As is common we assume this to have zero mean, so the multi-task GP is fully specified by the 
covariances (f T {x)f T ' (x 1 )) — C(t, x, t' , x'). For this covariance we take the flexible form from J5), 
(f T (x)f T i (x 1 )) = D TT iC(x, x'). Here C(x, x') determines the covariance between function values 
at different input points, encoding "spatial" behaviour such as smoothness and the lengthscale(s) 
over which the functions vary, while the matrix D is a free-form inter-task covariance matrix. 

One of the attractions of GPs for regression is that, even though they are non-parametric models 
with (in general) an infinite number of degrees of freedom, predictions can be made in closed form, 
see e.g. JT]. For a test point x for task r, one would predict as output the mean of f T (x) over 
the (Gaussian) posterior, which is y T K~ 1 k T (x). Here K is the n x n Gram matrix with entries 
Kim — D T(Tm C(xi, x m ) + a% 8i m , while k T (x) is a vector with the n entries k T j = D T(T C(xg, x). 
The error bar would be taken as the square root of the posterior variance of f T {x), which is 

V T {x) = D TT C{x, x) - kJ(x)K- 1 k T (x) (1) 

The learning curve for task r is defined as the mean-squared prediction error, averaged over the 
location of test input x and over all data sets with a specified number of examples for each task, say 
rt\ for task 1 and so on. As is standard in learning curve analysis we consider a matched scenario 
where the training outputs yg are generated from the same prior and noise model that we use for 
inference. In this case the mean-squared prediction error e T is the Bayes error, and is given by the 
average posterior variance [1], i.e. e T — (V T (x)) x . To obtain the learning curve this is averaged 
over the location of the training inputs xg. t T — (e T ). This average presents the main challenge 
for learning curve prediction because the training inputs feature in a highly nonlinear way in V T (x). 
Note that the training outputs, on the other hand, do not appear in the posterior variance V T (x) and 
so do not need to be averaged over. 

We now want to write the Bayes error e T in a form convenient for performing, at least approxi- 
mately, the averages required for the learning curve. Assume that all training inputs xg, and also the 
test input x, are drawn from the same distribution P(x). One can decompose the input-dependent 
part of the covariance function into eigenfunctions relative to P{x), according to C(x,x') = 
J2i ^i4>i( x )4>i( x ')- The eigenfunctions are defined by the condition (C(x, x')<pi(x')) x i — \i<f>i(x) 
and can be chosen to be orthonormal with respect to P(x), (<f>i(x)(f)j(x)) x — <Jy. The sum over i 
here is in general infinite (unless the covariance function is degenerate, as e.g. for the dot product 
kernel C(x, x') = x ■ x'). To make the algebra below as simple as possible, we let the eigenvalues 
Xi be arranged in decreasing order and truncate the sum to the finite range i — 1, . . . , M; M is then 
some large effective feature space dimension and can be taken to infinity at the end. 
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In terms of the above eigenfunction decomposition, the Gram matrix has elements 

i i,T,j,T r 

or in matrix form K = Sl/i\t' T + S where S is the diagonal matrix from the noise variances and 

^Lir = 5 TUT 4>i{xi), L iTt j T > = Xi6ijD TT > (2) 

Here ^ has its second index ranging over M (number of kernel eigenvalues) times T (number of 
tasks) values; L is a square matrix of this size. In Rronecker (tensor) product notation, L = D ® A 
if we define A as the diagonal matrix with entries Aj<5y. The Rronecker product is convenient for 
the simplifications below; we will use that for generic square matrices, (A <Ei B)(A' <Ei B') = 
(AA') ® (BB 1 ), (A ® B)^ 1 = A- 1 ®B- 1 , and tr (A ® B) = (tr A)(tr B). In thinking about 
the mathematical expressions, it is often easier to picture Rronecker products over feature spaces 
and tasks as block matrices. For example, L can then be viewed as consisting ofTxT blocks, each 
of which is proportional to A. 

To calculate the Bayes error, we need to average the posterior variance V T (x) over the test input x. 
The first term in ([T]i then becomes D TT (C(x, x)) — D TT tr A. In the second one, we need to average 

(C(xi, x)C(x, x m )) x D TmT 
= D TTl 2J \i\j(f>i{xt){(t)i(x)<j)j(x)) x (j)j{x m )D TmT 

ij 

= S ' D TT ''i l ii T i XiXjSij^mjr" D T >> T 
i,T',j,r" 

In matrix form this is {k T (x)kJ{x)) x = ^[{De T ejD) ® A 2 ]* T = \I*M T \I> T Here the last 
equality defines M T , and we have denoted by e T the T-dimensional vector with r-th component 
equal to one and all others zero. Multiplying by the inverse Gram matrix K 1 and taking the trace 
gives the average of the second term in ([T}; combining with the first gives the Bayes error on task r 

i T = (V T (x)) x = L> rr trA-tr*M r * T (* J L* T + S)" 1 

Applying the Woodbury identity and re-arranging yields 

e T = L> Tr trA-trM T * T S- 1 *(/ + i* T S- 1 *)- 1 
= L> TT trA-trM r i- 1 [/- (/ + J L* T S- 1 *)- 1 ] 

But 

trMrL- 1 = tr {[{De T eJ D) (g> A 2 }[D (g> A}- 1 } 

= tr{[De r e^] <x> A} = ejDe T tr A = D TT tr A 

so the first and second terms in the expression for e T cancel and one has 

e T = tr MrL- 1 (I + LV T -£- 1 '&)- 1 = trL^MrL-^L- 1 + ^E^ 1 *)" 1 
= tr [D ® A]- l [{De T e$D) <g> A 2 }[D <x> A]" 1 ^" 1 + # T 5r x *) -1 
= tr[e T e^(g) ^(X" 1 + ^S" 1 *)" 1 

The matrix in square brackets in the last line is just a projector P T onto task t; thought of as a matrix 
of T x T blocks (each of size M x M), this has an identity matrix in the (r, r) block while all other 
blocks are zero. We can therefore write, finally, for the Bayes error on task r, 

e T = tiP T (L~ 1 + * T S- 1 *)- 1 (3) 

Because X is diagonal and given the definition (|2ji of \&, the matrix \I> T £ _1 \l r is a sum of contri- 
butions from the individual training examples I = 1, . . . , n. This will be important for deriving the 
learning curve approximation below. We note in passing that, because ^ T P T = I, the sum of the 
Bayes errors on all tasks is J2 T i T — tr (L^ 1 + \I> T £ _1 \(') _1 , in close analogy to the corresponding 
expression for the single-task case lfl3l . 
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3 Learning curve prediction 



To obtain the learning curve e T = (e T ), we now need to carry out the average (...) over the training 
inputs. To help with this, we can extend an approach for the single-task scenario |fl3l and define a 
response or resolvent matrix Q = (L^ 1 + \I> T E _1 \I> + ^ T ti T P T ) _1 with auxiliary parameters v T 
that will be set back to zero at the end. One can then ask how G = (Q) and hence e r < = tr P T >G 
changes with the number n T of training points for task r. Adding an example at position x for task 
r increases \I' T X -ll 4' by <j~ 2 <j) T 4>J , where </> T has elements (</> T ); r < = <j>i(x)S TT /. Evaluating the 
difference (Q^ 1 + a^^rcjij)^ 1 — Q with the help of the Woodbury identity and approximating it 
with a derivative gives 

This needs to be averaged over the new example and all previous ones. If we approximate by 
averaging numerator and denominator separately we get 



dG 1 dG 



dn T a 2 + tr P T G dv T 



(4) 



Here we have exploited for the average over x that the matrix (4> T 4>^) X has («,t'), (j, r")-entry 
(4>i(x)cf>j(x)) x 5 TT 'd TT " = Sij5 TT rS TT ", hence simply (4> T 4>J) X = P T . We have also used the 
auxiliary parameters to rewrite — (QP T Q) = d(Q)/dv T = dG/dv T . Finally, multiplying Q by P T i 
and taking the trace gives the set of quasi-linear partial differential equations 

de T , _ 1 de T , 
dn T a 2 + e T dv T 

The remaining task is now to find the functions e T (n l5 . . . , nx, «i, • • • , vt) by solving these differ- 
ential equations. We initially attempted to do this by tracking the e T as examples are added one 
task at a time, but the derivation is laborious already for T = 2 and becomes prohibitive beyond. 
Far more elegant is to adapt the method of characteristics to the present case. We need to find a 
2T-dimensional surface in the 3T-dimensional space (m, . . . , ny, v±, . . . ,Vt, eij • • • , £t), which is 
specified by the T functions e r (. . .). A small change [5n\, . . . , Snr, Svi, . . . , Svt, 6e±, . . . , Sex) in 
all 3T coordinates is tangential to this surface if it obeys the T constraints (one for each r') 

Se T , = S~] (^-Srir + ^-Sv T 

T 

From (|5j, one sees that this condition is satisfied whenever Se T = and Sn T = —8v T {<7 2 + e T ) 
It follows that all the characteristic curves given by e T (t) = e T = const., v T (t) = w Ti o(l — t), 
n T (t) = v T fl(a 2 + e T .o) t for t G [0, 1] are tangential to the solution surface for all t, so lie within 
this surface if the initial point at t = does. Because at t = there are no training examples 
(n T (0) = 0), this initial condition is satisfied by setting 

e T , = trP r (^L- 1 + «r',o-Pr'J 

Because e T (t) is constant along the characteristic curve, we get by equating the values at t = and 
t = 1 

e T ,o = trP T r 1 + ^ v Tlfi P T ^j = e T {{n T < = v T >. a {a 2 T , + e T ',o)}, {v T ' = 0}) 
Expressing v T > in terms of n T ' gives then 

e T = tr P T (l-i + -T^—Pr) (6) 



This is our main result: a closed set of T self-consistency equations for the average Bayes errors 
e T . Given L as defined by the eigenvalues A, of the covariance function, the noise levels a 2 , and the 
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number of examples n T for each task, it is straightforward to solve these equations numerically to 
find the average Bayes error e T for each task. 

The r.h.s. of |6]i is easiest to evaluate if we view the matrix inside the brackets as consisting of 
M x M blocks of size T x T (which is the reverse of the picture we have used so far). The matrix is 
then block diagonal, with the blocks corresponding to different eigenvalues A^. Explicitly, because 
L^ 1 = D 1 ® A -1 , one has 

e T =Y,(K l D-'+^g{{-^- })) 1 (7) 



4 Results and discussion 



We now consider the consequences of the approximate prediction (j7]i for multi-task learning curves 
in GP regression. A trivial special case is the one of uncorrected tasks, where D is diagonal. Here 
one recovers T separate equations for the individual tasks as expected, which have the same form as 
for single-task learning lfl3l . 



4.1 Pure transfer learning 

Consider now the case of pure transfer learning, where one is learning a task of interest (say r = 1) 
purely from examples for other tasks. What is the lowest average Bayes error that can be obtained? 
Somewhat more generally, suppose we have no examples for the first Tq tasks, nx = • • • = nx a = 0, 
but a large number of examples for the remaining T\ = T — To tasks. Denote E = D~ Y and write 
this in block form as 




Now multiply by A^ 1 and add in the lower right block a diagonal matrix N = dmg{{n T /{a 2 + 
e T )} r= T +i ,....t)- The matrix inverse in ^ then has top left block X^Eqq + Eqq Eox(\iN + 
En — EqxEqq Eqx)~ 1 Eq 1 Eq ]. As the number of examples for the last T\ tasks grows, so do all 
(diagonal) elements of AT. In the limit only the term \iE^ survives, and summing over i gives 
6i = tr A(^ 00 1 )ii = (C(x, x)) (Eqq)xx- The Bayes error on task 1 cannot become lower than this, 

P.cing a limit on the benefits of pure transfer learning. That this prediction of the approximation 
for such a lower limit is correct can also be checked directly: once the last Xi tasks f T {x) 
= To + 1, . . . T) have been learn perfectly, the posterior over the first To functions is, by standard 
Gaussian conditioning, aGP with covariance C(x,x'){Eqq)~ x . Averaging the posterior variance of 
fi (x) then gives the Bayes error on task 1 as e\ — (C (x,x)) {Eqq ) h, as found earlier. 

This analysis can be extended to the case where there are some examples available also for the first 
T tasks. One finds for the generalization errors on these tasks the prediction ^ with D 1 replaced 
by £^oo- This is again in line with the above form of the GP posterior after perfect learning of the 
remaining T\ tasks. 



4.2 Two tasks 



We next analyse how well the approxiation (f7| does in predicting multi-task learning curves for 
T = 2 tasks. Here we have the work of Chai [21] as a baseline, and as there we choose 



D = 



1 P 
P 1 



The diagonal elements are fixed to unity, as in a practical application where one would scale both 
task functions fx(%) and f2(x) to unit variance; the degree of correlation of the tasks is controlled 
by p. We fix tt2 = U2 /n and plot learning curves against n. In numerical simulations we ensure 
integer values of nx and n-2 by setting n-2 = L n7r 2j, n\ = n — ri2\ for evaluation of (fvj we use 
directly n% = mr2, tlx = n(l — 7ra). For simplicity we consider equal noise levels Oi=o\ = a 2 . 

As regards the covariance function and input distribution, we analyse first the scenario studied 
in ETI : a squared exponential (SE) kernel C(x, x') = exp[— (x — x') 2 /(2l 2 )] with lengthscale 
I, and one-dimensional inputs x with a Gaussian distribution JV(0, 1/12). The kernel eigenvalues A; 
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n n n 



Figure 1 : Average Bayes error for task 1 for two-task GP regression with kernel lengthscale / = 0.01, 
noise level a 2 = 0.05 and a fraction tt2 = 0.75 of examples for task 2. Solid lines: numerical 
simulations; dashed lines: approximation f7j. Task correlation p 2 = 0, 0.25, 0.5, 0.75, 1 from top to 
bottom. Left: SE covariance function, Gaussian input distribution. Middle: SE covariance, uniform 
inputs. Right: OU covariance, uniform inputs. Log-log plots (insets) show tendency of asymptotic 
uselessness, i.e. bunching of the p < 1 curves towards the one for p = 0; this effect is strongest for 
learning of smooth functions (left and middle). 

are known explicitly from [22 J and decay exponentially with i. Figure [TJleft) compares numerically 
simulated learning curves with the predictions for €\, the average Bayes error on task 1, from Q. 
Five pairs of curves are shown, for p 2 — 0, 0.25, 0.5, 0.75, 1. Note that the two extreme values 
represent single-task limits, where examples from task 2 are either ignored (p = 0) or effectively 
treated as being from task 1 (p = 1). Our predictions lie generally below the true learning curves, 
but qualitatively represent the trends well, in particular the variation with p 2 . The curves for the dif- 
ferent p 2 values are fairly evenly spaced vertically for small number of examples, n, corresponding 
to a linear dependence on p 2 . As n increases, however, the learning curves for p < 1 start to bunch 
together and separate from the one for the fully correlated case (p — 1). The approximation (j7]) 
correctly captures this behaviour, which is discussed in more detail below. 

Figure [TJ middle) has analogous results for the case of inputs x uniformly distributed on the interval 
[0, 1]; the Xi here decay exponentially with i 2 ifTTl . Quantitative agreement between simulations 
and predictions is better for this case. The discussion in [17| suggests that this is because the 
approximation method we have used implicitly neglects spatial variation of the dataset-averaged 
posterior variance (V T (x)}; but for a uniform input distribution this variation will be weak except 
near the ends of the input range [0, 1]. Figure mright) displays similar results for an OU kernel 
C(x,x') = exp(— |a; — x'\/l), showing that our predictions also work well when learning rough 
(nowhere differentiable) functions. 

4.3 Asymptotic uselessness 

The two-task results above suggest that multi-task learning is less useful asymptotically: when the 
number of training examples n is large, the learning curves seem to bunch towards the curve for 
p = 0, where task 2 examples are ignored, except when the two tasks are fully correlated (p = 1). 
We now study this effect. 

When the number of examples for all tasks becomes large, the Bayes errors e T will become small 
and eventually be negligible compared to the noise variances a 2 in One then has an explicit 
prediction for each e r , without solving T self-consistency equations. If we write, for T tasks, 
n T = nir T with ir T the fraction of examples for task r, and set 7 T = ir T /a 2 , then for large n 

= E ( (K'd- 1 + nr)" 1 = E^r-^A^rVW/ 2 )- 1 + m\-^' 2 ) TT (8) 

where T = diag(7i , . . . , 7r). Using an eigendecomposition of the symmetric matrix r 1 / 2 £>r 1 / 2 = 
Ea=i $aV a Va, one then shows in a few lines that (|8| can be written as 

Er^T^Eatvf^H) (9) 
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Figure 2: Left: Bayes error (parameters as in Fig. [gleft), with n = 500) vs p 2 . To focus on the 
error reduction with p, r = [ei(p) — ei(l)]/[ei(0) — ei(l)] is shown. Circles: simulations; solid 
line: predictions from (jTjl. Other lines: predictions for larger n, showing the approach to asymptotic 
uselessness in multi-task learning of smooth functions. Inset: Analogous results for rough functions 
(parameters as in FigJiTright)). Right: Learning curve for many-task learning (T — 200, parameters 
otherwise as in Fig.luTeft) except p 2 = 0.8). Notice the bend around ei = 1 — p = 0.106. Solid 
line: simulations (steps arise because we chose to allocate examples to tasks in order t = 1, . . . , T 
rather than randomly); dashed line: predictions from Q. Inset: Predictions for T = 1000, with 
asymptotic forms e = 1 — p + pi and e = (1 — p)e for the two learning stages shown as solid lines. 



where g(h) — tr(A _1 + h)^ 1 = ^ i (A i _1 + h)^ 1 and v a . r is the r-th component of the a-th 
eigenvector v a . This is the general asymptotic form of our prediction for the average Bayes error 
for task r. 

To get a more explicit result, consider the case where sample functions from the GP prior have 
(mean-square) derivatives up to order r. The kernel eigenvalues \i then decay a^]i~( 2r+2 ) for large 
i, and using arguments from ifTTll one deduces that g(h) ~ h~ a for large h, with a = (2?' + l)/(2r + 
2). In ^ we can then write, for large n, g(nS a ) (S a /^f T )^ a g(n r y T ) and hence 

e T a 5(«7r){EaK,.) 2 (^/7r) 1 - Q } (10) 

When there is only a single task, Si — 71 and this expression reduces to ei = g(nji) = g{ni/a 2 ). 
Thus g(nj T ) = g(n T /cr 2 ) is the error we would get by ignoring all examples from tasks other than r, 



and the term in { . . .} in ( 10 1 gives the "multi-task gain", i.e. the, factor by which the error is reduced 
because of examples from other tasks. (The absolute error reduction always vanishes trivially for 
n — » 00, along with the errors themselves.) One observation can be made directly. Learning of very 
smooth functions, as defined e.g. by the SE kernel, corresponds to r — > 00 and hence a — > 1, so 
the multi-task gain tends to unity: multi-task learning is asymptotically useless. The only exception 
occurs when some of the tasks are fully correlated, because one or more of the eigenvalues 5 a of 
r 1 /2£ )r i/2 will then be zero 

left) shows this effect in action, plotting Bayes error against p 2 for the two-task setting of 
left) with n = 500. Our predictions capture the nonlinear dependence on p 2 quite well, 
though the effect is somewhat weaker in the simulations. For larger n the predictions approach a 
curve that is constant for p < 1, signifying negligible improvement from multi-task learning except 
at p = 1. It is worth contrasting this with the lower bound from [21 1, which is linear in p 2 . While 
this provides a very good approximation to the learning curves for moderate n [21], our results here 
show that asymptotically this bound can become very loose. 

When predicting rough functions, there is some asymptotic improvement to be had from multi-task 
learning, though again the multi-task gain is nonlinear in p 2 : see Fig. |2jleft, inset) for the OU case, 
which has r = 1). A simple expression for the gain can be obtained in the limit of many tasks, to 
which we turn next. 




1 See the discussion of Sacks- Ylvisaker conditions in e.g. 1 1 1; we consider one-dimensional inputs here 
though the discussion can be generalized. 
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4.4 Many tasks 



We assume as for the two-task case that all inter-task correlations, D T T i with r 7^ r', are equal to 
p, while D T/T = 1. This setup was used e.g. in 11231 . and can be interpreted as each task having a 
component proportional to ^fp of a shared latent function, with an independent task-specific signal 
in addition. We assume for simplicity that we have the same number n T — n/T of examples for 
each task, and that all noise levels are the same, a 2 — a 2 . Then also all Bayes errors e T = e will be 
the same. Carrying out the matrix inverses in (|7]i explicitly, one can then write this equation as 

e = g T (n/(a 2 + e),p) (11) 
where gr{h, p) is related to the single-task function g(h) from above by 

g T (h, P ) = ^(1 - p)g(h(l - p)/T) + (p + g(h[p + (1 - p)/T]) (12) 

Now consider the limit T —> 00 of many tasks. If n and hence h = n/{a 2 + e) is kept 
fixed, grp{h,p) — > (1 — p) + pg(hp); here we have taken g(0) = 1 which corresponds to 



tr A = (C(x, x)) x = 1 as in the examples above. One can then deduce from (Hi that the Bayes 
error for any task will have the form e = (1 — p) +pe, where 1 decays from one to zero with increas- 
ing n as for a single task, but with an effective noise level a 2 = (1 — p + a 2 ) / p. Remarkably, then, 
even though here n/T -> so that for most tasks no examples have been seen, the Bayes error for 
each task decreases by "collective learning" to a plateau of height 1 — p. The remaining decay of e to 
zero happens only once n becomes of order T. Here one can show, by taking T — >• 00 at fixed h/T 



12 1 and inserting into (111, that e = (1 — p)e where e again decays as for a single task but with 
an effective number of examples n = n/T and effective noise level ct 2 /(1 — p). This final stage of 
learning therefore happens only when each task has seen a considerable number of exampes n/T. 
Fig. fright) validates these predictions against simulations, for a number of tasks (T = 200) that 
is in the same ballpark as in the many-tasks application example of [24] . The inset for T = 1000 
shows clearly how the two learning curve stages separate as T becomes larger. 

Finally we can come back to the multi-task gain in the asymptotic stage of learning. For GP priors 
with sample functions with derivatives up to order r as before, the function e from above will decay 
as {n/cr 2 )~ a ; since e = (1 — p)e and a 2 — cr 2 /(l — p), the Bayes error e is then proportional 
to (1 — p) 1_Q - This multi-task gain again approaches unity for p < 1 for smooth functions (a = 
(2r + 1)/ (2r + 2) — > 1). Interestingly, for rough functions (a < 1), the multi-task gain decreases 
for small p 2 as 1 — (1 — a)^/]} 2 and so always lies below a linear dependence on p 2 initially. This 
shows that a linear-in-p 2 lower error bound cannot generally apply to T > 2 tasks, and indeed one 
can verify that the derivation in [21 ] does not extend to this case. 



5 Conclusion 

We have derived an approximate prediction for learning curves in multi-task GP regression, valid 
for arbitrary inter-task correlation matrices D. This can be evaluated explicitly knowing only the 
kernel eigenvalues, without sampling or recourse to single-task learning curves. The approximation 
shows that pure transfer learning has a simple lower error bound, and provides a good qualitative 
account of numerically simulated learning curves. Because it can be used to study the asymptotic 
behaviour for large training sets, it allowed us to show that multi-task learning can become asymp- 
totically useless: when learning smooth functions it reduces the asymptotic Bayes error only if tasks 
are fully correlated. For the limit of many tasks we found that, remarkably, some initial "collective 
learning" is possible even when most tasks have not seen examples. A much slower second learning 
stage then requires many examples per task. The asymptotic regime of this also showed explicitly 
that a lower error bound that is linear in p 2 , the square of the inter-task correlation, is applicable 
only to the two-task setting T = 2. 

In future work it would be interesting to use our general result to investigate in more detail the con- 
sequences of specific choices for the inter-task correlations D, e.g. to represent a lower-dimensional 
latent factor structure. One could also try to deploy similar approximation methods to study the case 
of model mismatch, where the inter-task correlations D would have to be learned from data. More 
challenging, but worthwhile, would be an extension to multi-task covariance functions where task 
and input-space correlations to not factorize. 
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